SU{4:) Fermi Liquid State and Spin Filtering in a Double Quantum Dot System 
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We study a symmetrical double quantum dot (DD) system with strong capacitive inter-dot coupling 
using renormalization group methods. The dots are attached to separate leads, and there can be a 
weak tunneling between them. In the regime where there is a single electron on the DD the low- 
energy behavior is characterized by an S'f/(4)-symmetric Fermi liquid theory with entangled spin 
and charge Kondo correlations and a phase shift tt/A. Application of an external magnetic field gives 
rise to a large magneto-conductance and a crossover to a purely charge Kondo state in the charge 
sector with 5*17(2) symmetry. In a four lead setup we find perfectly spin polarized transmission. 



PACS numbers; 75.20.Hr, 71.27.-ha, 72.15.Qm 

Introduction. — Quantum dots are one of the most ba- 
sic building blocks of mesoscopic circuits [1]. In many 
respects quantum dots act as large complex atoms cou- 
pled to conducting leads that are used to study transport. 
The physical properties of these dots depend essentially 
on the level spacing and precise form of the coupling to 
the leads: They can exhibit Coulomb blockade phenom- 
ena [2] , build up correlated Kondo- like states of various 
kinds [3-5], or develop conductance fluctuations. 

The simplest mesoscopic circuits that go beyond sin- 
gle dot devices in their complexity are double dot (DD) 
devices (see Fig. 1). These 'artificial molecules' have 
been extensively studied both theoretically [6-11] and 
experimentally [12-15]: They may give rise to stochastic 
Coulomb blockade [6] and peak splitting [7,12], can be 
used as single electron pumps [1], were proposed to mea- 
sure high frequency quantum noise [11], and are building 
blocks for more complicated mesoscopic devices such as 
turn-stiles or cellular automata [16]. DD's also have in- 
teresting degeneracy points where quantum fluctuations 
may lead to unusual strongly correlated states [17]. 

In the present paper we focus our attention to small 
semiconducting DDs with large inter-dot capacitance 
[10,17]. We consider the regime where the gate volt- 
ages V± are such that the lowest lying charging states, 
(7i4_,7i_) — (0,1) and (1,0), are almost degenerate: 
i?(l, 0) — E{0, 1) « [n± = # of extra electrons on dot 
'±', and E[n-^^, n_) is measured from the common chem- 
ical potential of the two leads]. We consider the simplest, 
most common case where the states (1, 0) and (0, 1) have 
both spin S — 1/2, associated with the extra electron on 
the dots. Then at energies below the charging energy of 
the DD, Ec = min{i;(l, 1) - E{Q, 1), S(0, 0) - E{Q, 1)}, 
the dynamics of the DD is restricted to the subspace 
{S'- =±1/2; n+-n_ =±1}. 

As we discuss below, quantum fluctuations between 
these four quantum states of the DD generate an un- 
usual strongly correlated Fermi liquid state, where the 
spin and charge degrees of freedom of the DD are totally 



entangled. We show that this state possesses an SU{A) 
symmetry corresponding to the total internal degrees of 
freedom of the DD, and is characterized by a phase shift 
(5 = 7r/4. This phase shift can be measured by integrating 
the DD device in an Aharonov-Bohm interferometer [18]. 
Application of an external field on the DD suppresses spin 
fluctuations. However, charge fluctuations are unaffected 
by the magnetic field and still give rise to a Kondo effect 
in the charge (orbital) sector [10,17,19]. We show that 
in a four lead setup this latter state gives rise to an al- 
most totally spin polarized current through the DD with 
a field-independent conductance G « /h. The conduc- 
tance across the dots, on the other hand, shows a large 
negative magneto-resistance at T = temperature. 

Model. — Let us first discuss the arrangement in Fig. 1 . 
At energies below Ec we describe the isolated DD in 
terms of the orbital pseudospin = (n_|_ — n^)/2 = i^: 



(1) 



Hdot = ~SE T"" -tT"" - B 



The term proportional to describes the energy differ- 
ence of the two charge states [5E = E{1,0) — £'(0, 1) ~ 
V+ — for a fully symmetrical system], while t <^ Ec 
is the tunneling amplitude between them. The last term 
stands for the Zeeman-splitting due to an applied local 
magnetic field in the z direction. We are interested in the 
regime, where — despite the large capacitive coupling, — 
the tunneling between the dots is small. Furthermore, 
one needs a large enough single particle level spacing A 
on the dots. Both conditions can be satisfied by making 
small dots [20], which are close together or capacitively 
coupled to a common top-gate electrode [21]. 
The leads are described by the Hamiltonian: 



|e|<-D 



J2 

|£|<D 



(2) 



where al^^ (Oso--) creates an electron in the right (left) 
lead with energy e and spin a, D ^ Tmii{Ec, A} = 1 is 



a cut-off, and {a| 



'r'} = drr'Saa'Sie - e'). 
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FIG. 1. Top: Schematics of the DD device. Bottom: Vir- 
tual process leading to 'spin-flip assisted tunneling' as de- 
scribed in Eq. (4) 

To determine the effective DD lead coupling we have 
to consider virtual charge fluctuations to the excited 
states with n+ + ri- =0 and 2, generated by tunneling 
from the leads to the dots. By second order perturbation 
theory in the lead-dot tunneling we obtain the following 
effective Hamiltonian: 

i^Kondo = ^J+P+S{iP^ap+iP) + ^J-P-S{i;^ap-tP), (3) 

i?a«sist = Q± [T+S{tP^aT-tP) + h.c.] , (4) 

Horb = i (l4T"(VV"V) + V1[T+(VV-V) + h.c.]). (5) 

where Vctt = J de agar, and a and r denote the spin 
and orbital pseudospin of the electrons (cr =t,i; t = 
= ±1). The operators P± = (1 ± 2r^)/2 and p± = 
(l±T^)/2 project out the DD states (1, 0) and (0, 1), and 
the right/left lead channels, respectively. 

In the limit of small dot-lead tunneling the dimension- 
less exchange couplings are J± ~ T±/Ec with r± the 
tunneling rate to the right (left) lead [22]. The 'spin- flip 
assisted tunneling' Qx ~ -\/r+r3/^c in Eq. (4) gives 
simultaneous spin- and pseudospin-flip scattering and is 
produced by virtual processes depicted in the lower part 
of Fig. 1, while the spin- independent parts of such virtual 
processes lead to the orbital Kondo term in Eq. (5) with 
similar amplitudes. 

We first focus on the case of a fully symmetrical DD. 
Then the sum of Eqs. (3) and (4) can be rewritten as: 

-ffKondo + -ffassist = S{^'^ Stjj) (6) 

+ Q.T'Sii^^T'ai;) + Qi_ (r+5(V'V-aV) + h.c.) , 

where J = Qz = {J+ + J-)/4. The coupHngs in Eqs. (3- 
5) are not entirely independent, but are related by the 
constraints Vj_ = Q± and J = Qz- 



Scaling Analysis. — The perturbative scaling analysis 
follows that of a related model in Ref. [23]. In the per- 
turbative RG one performs the scaling by integrating out 
conduction electrons with energy larger than a scale scale 
D <^ D, and thus obtains an effective Hamiltonian that 
describes the physics at energies D. For SE = t = B = 0, 
in the leading logarithmic approximation we find that all 
couplings diverge at the Kondo temperature where 
the perturbative scaling breaks down. Nevertheless, the 
structure of the divergent couplings suggests that at low 
energies J = V± = Vz = Q± = Qz- Thus at small ener- 
gies, — apart from a trivial potential scattering — the 
effective model is a remarkably simple SU(4) symmetrical 
exchange model: 

Heff(r-o) = j J2 V'iV'/3 , (7) 

a,/3=l,..,4 

where the index a labels the four combinations of possible 
spin and pseudospin indices, and the |q:)'s denote the four 

states of the DD. This indeed can be more rigorously 
proven using strong coupling expansion, conformal field 
theory, and large / (flavor) expansion techniques [24-26], 
and is also confirmed by out numerical computations. 

Numerical Renormalization Group (NRG). — To ac- 
cess the low-energy physics of the DD, we used Wilson's 
NRG approach [27] . In this method one defines a series of 
rescaled Hamiltonians, Hn, related by the relation [27]: 

iJjV+l = A^/^i/jv + {fl,arfN+l,ar + h.c) , (8) 

(TT 

where foar ^ tpar/V^ and Hq = 2Ai/2/(l + A) Hint with 
A ~ 3 as discretization parameter, and ^at « 1. (For the 
definition of see Ref. [27].) We have defined Hint = 
Hdot + -ffKondo + -^assist + i^orb . The Original Hamilto- 
nian is related to the Hn's as H = hmiv^oo <jJnHn with 
LUN = A"(^+i)/2(l + A)/2. Solving Eq. (8) iteratively we 
can then use the eigenstates of Hn to calculate physical 
quantities at a scale T,ui u)n- 

Results. — First let us consider the case if^ot = 0. 

Fixed point structure- — The finite size spectrum pro- 
duced by the NRG procedure contains a lot of informa- 
tion. Among others, we can identify the structure of the 
low-energy effective Hamiltonian from it [27], and also 
determine all scattering phase shifts. 

In particular, we find that for SE = t — B — the 
entire finite size spectrum can be understood as a sum 
of four independent, spinless chiral fermion spectra with 
phase shifts 6 = 7r/4. This phase shift is characteristic 
for the SU{4) Hamiltonian, Eq. (7), and simply follows 
from the Friedel sum rule [24]. Application of an exter- 
nal magnetic field B to the DD gradually shifts 6 to the 
values (5| ^ 7r/2 and ^ [28]. 

Spectral functions. — To learn more about the dynam- 
ics of the DD we computed at T = the spin spec- 
tral function Qg{u!) = — (l/7r)Im{x|(a;)}, and pseudospin 
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FIG. 2. T — spin and pseudospin spectral functions for 
J — Qz = Vz = 0.14, V± — Q± — 0.13, and various val- 
ues of B EE B/T^\ For B = both the spin and pseudospin 
spectral functions exhibit ^ to behavior below the Kondo tem- 
10""^. Applying a magnetic field the situa- 



perature T^-* 



tion changes: The B > magnetic field destroys the spin 
Kondo correlations and leads to a purely orbital Kondo-effect. 

spectral function q'~1-{lo) = — (l/7r)Im{x^(cj)} by the den- 
sity matrix NRG method [29]. 

At -B = the various spectral functions exhibit a peak 
at the same energy, , corresponding to the formation 

of the SU{A) symmetric state (see Fig. 2). Below T^-* 
all spectral functions become Hnear, characteristic to a 
Fermi liquid state with local spin and pseudospin sus- 
ceptibilities Xs ~ At ~ 1 /Tk'^ , where the SU (4) " hyper- 
spin" of the dot electron (formed by {t +,i +,t ^li 
components) is screened by the lead electrons. 

Now let us consider the case i?dot 7^ 0. In a large mag- 
netic field, T^^ <C B, spin-flip processes are suppressed: 
The spin spectral function therefore shows only a Schot- 
tky anomaly dX lo ^ B. Nevertheless, the couplings VI 
and Vz still generate a purely orbital Kondo state in the 
spin channel with the same orientation as the DD spin, 
with a reduced Kondo temperature Tx{B) < and 
a corresponding phase shift 5^ ~ 7r/2. 

Due to the spin-pseudospin symmetric structure of the 
Hamiltonian, Eq. (6), the opposite effect occurs for a 
large 6E: In that limit the charge is localized on one side 
of the DD, charge fluctuations are suppressed, and the 
system scales to a spin Kondo problem. A large tun- 
neling, t > T^'' is also expected to lead to a somewhat 
similar effect, though the conductance through the DD 
behaves very differently in the two cases [28] . 

DC Conductivity. — First we focus on the conductivity 
across the DD assuming a small tunneling t. Then we 
can assume that the two dots are in equilibrium with 
the leads connected to them, and we can compute the 
induced current perturbatively in t. A simple calculation 
yields the following formula: [30] 
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FIG. 3. The T = conductance of DD system at a; = 
for J = Qz = 0.14, Vj_ = Q_L = 0.13, = 0.14 and different 
magnetic field values. Inset: small B limit of the conductance. 



G = 



r lim 



(9) 



The normalized DC conductance at T = temperature 
is shown in Fig 3. Below the orbital Kondo temperature 
gy{uj) ^ uj/T^{B), leading to a dimensionless conduc- 
tance ~ [t/TK{B)]^ . However, Tk{B) strongly decreases 
with increasing B implying a large negative magneto- 
resistance in the T = DC conductance. This effect is re- 
lated to the correlation between spin and orbital degrees 
of freedom. We have to emphasize that the simple consid- 
erations above only apply in the regime t <C Tk{B). For 
larger values of t a more complete calculation is needed. 

Having extracted the phase shifts from the NRG spec- 
tra, we can construct the scattering matrix in more gen- 
eral geometries too and compute the T = conductance 
using the Landauer-Buttiker formula [28,31,32]. In the 
perfectly symmetrical two terminal four lead setup of 
Fig. 4 with 5E — t — Q, e.g., the DC conductance is 
Gi3 = \GQ{sin^{5i{B)) -f sin^{5^{B))), where Gq = 
2e^ jh is the quantum conductance. By the Friedel sum 
rule b^{B) =tt/2-5i{B), and thus Gi3(r = 0) = Gg/2, 
independently of B. However, the polarization of the 
transmitted current, P = 2sin^((5|) — 1 tends rapidly to 
1 as i? > rj?"* , and the DD thereby acts as a perfect spin 

filter at r = with B > , and could also serve as 
a spin pump. For a typical Tk ~ 0.5K and a g-factor 
g Ri —0.4 as in GaAs, e.g. a field of 2.5T would give 
a 97% polarized current. This efficiency is comparable 
to other spin filter designs [33]. By lowering Tk, much 
higher polarizations could be obtained. 

Robustness. — Since the spin S" and pseudospin 
are both marginal operators at the S'C/(4) fixed point 
[25], we conclude that the S'C/(4) behavior is stable in 
the sense that a small but finite value of SE, B, t <^ 
T^' will lead only to small changes in physical properties 
like the phase shifts. The anisotropy of the couplings is 
also irrelevant in the RG sense [25,26], and the role of 
J_ 7^ J4- symmetry breaking is only to renormalize the 
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FIG. 4. Top: Magnetic field dependence of the phase shifts 
for t = SE « 0. Bottom: Corresponding T = dimensionless 
conductance and spin polarization of the current in the four 
lead setup shown in the Top inset. 

bare value of dE^ which is a marginal perturbation itself. 
Therefore the S'(7(4) Fermi liquid state is robust under 
the conditions discussed in the Introduction. 

Experimental accessibility. — For our scenario it is cru- 
cial to have large enough charging energy and level spac- 
ing £'c,A > T^°^ > t. With today's technology it is 
possible to reach A 2 — 3K. The dot-dot capaci- 
tance (and thus Ec [7]) can be increased by changing 
the shape of the gate electrode separating the dots, us- 
ing a columnar geometry as in Refs. [19,34], where the 
two-dimensional dots are placed on the top of each other, 
or placing an additional electrode on the top of the DD 
device [21]. We could not find a closed expression for 
in the general case. However, for a symmetrical 
BB J ^ ^ V, ^ ^ Q, ^ T/27tEc, pro- 

vided that fluctuations to the (0, 0) state give the domi- 
nant contribution. Then we obtain T^-* ft! De~^^^"' and 

TxiB = oo) ft est. [rl^^f/D. Thus the value of J and 
thus can be tuned experimentally to a value similar 
to the single dot experiments. Indeed, an orbital Kondo 
effect has recently been observed [19]. 

Summary. — We have studied a DD system with large 
capacitive coupling close to its degeneracy point, in the 
Kondo regime. Using both scaling arguments and a non- 
perturbative NRG analysis, we showed that the simulta- 
neous appearance of the Kondo effect in the spin and 
charge sectors results in an S'C/(4) Fermi liquid ground 
state with a phase shift 7r/4. Upon applying an exter- 
nal magnetic field, the system crosses over to a purely 



charge Kondo state with a lower Tk- In a four-terminal 
setup, the DD could thus be used as a spin filter with 
high transmittance. We further predict a large serial 
magneto-conductance at T = 0. The SU{4:) behavior 
in this system is robust, and is experimentally accessible. 
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